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ABSTRACT 


Fracture  toughness  values  determined  using  shallow  cracked  single  edge  notch  bend, 
SB(B),  specimens  of  structural  thickness  are  useful  for  structural  integrity  assessments. 
However,  testing  standards  have  not  yet  incorporated  formulas  that  permit  evaluation  of 
J  and  CTOD  for  shallow  cracks  from  experimentally  measured  quantities  (i.e.  load,  crack 
mouth  opening  displacement  (CMOD),  and  load  line  displacement  (LLD)).  Pesults  from  two 
dimensional  plane  strain  finite-element  analyses  are  used  to  develop  J  and  CTOD  estima¬ 
tion  strategies  appropriate  for  application  to  both  shallow  and  deep  crack  SE(B)  specimens, 
Crack  depth  to  specimen  width  (a/W)  ratios  between  0.05  and  0.70  are  modelled  using 
Ramberg-Osgood  strain  hardening  exponents  (n)  between  4  and  50.  The  estimation  formu¬ 
las  divide  J  and  CTOD  into  small  scale  yielding(<SSy)  and  large  scale  yielding  (LSY)  compo¬ 
nents.  For  each  case,  the  SSY componer  t  is  determined  by  the  linear  elastic  stress  intensity 
factor,  Kj.  The  formulas  differ  in  evaluation  of  the  LSY  component.  The  techniques  consid¬ 
ered  include:  estimating  J  or  CTOD  from  plastic  work  based  on  load  line  displacement 
{Api  I  r  rn).  from  plastic  work  based  on  crack  mouth  opening  displacement  (A^;  I  cmoj^» 
from  tlmplastic  component  of  crack  mouth  opening  displacementiCMODjp/).  A.i  I  ca/odP^^* 
vides  the  most  accurate  J  estimation  possible.  Hie  finite-element  results  for  ali  conations 
investigated  fall  within  9%  of  the  following  formula: 

J  m  :  where  rjj_^  m  3.786  -  8.101  |p>  +  2.018|^j 

The  insensitivity  of  r}j_Q  to  strain  hardening  permits  J  estimation  for  any  material  with 
equal  accuracy.  Further,  estimating  J  from  CMOD  rather  than  LLD  eliminates  the  need 
to  measure  LLD,  thus  simplifying  the  test  procedure.  Alternate,  work  based  estimates  for 
J  and  CTOD  have  equivalent  accuracy  to  t^s  formula;  however  the  rj  coefTicients  in  these 
equations  depend  on  the  strain  hardening  coefficient.  CTOD  estimates  based  on  scalar  pro¬ 
portionality  of  CTOD/iry  and  CMODp/  are  highly  inaccurate,  especially  for  materials  with 
considerable  strain  hardening,  where  errors  up  to  38%  occur. 
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1.  INTRODUCTION 

Standardized  procedures  for  fracture  toughness  testing  require  both  sufficient  specimen  thickness  to 
insure  predominantly  plane  strain  conditions  at  the  crack  tip  and  a  crack  depth  of  at  least  half  the  spec¬ 
imen  width  [l-B].  Within  certain  limits  on  load  level  and  crack  growth,  these  restrictions  insure  the 
existence  of  veiy  severe  conditions  for  fracture  as  described  by  the  Hutchinson  Rice  Rosengren 
(HRR)  crack-tip  fields  [4,5].  These  conditions  make  the  applied  driving  force  needed  to  initiate  frac¬ 
ture  in  a  laboratory  specimen  lower  than  the  value  needed  to  initiate  fracture  in  common  civil  and 
marine  structures  where  such  severe  geometric  conditions  are  not  present.  As  a  consequence,  struc¬ 
tures  often  cany  greater  loads  without  failure  than  predicted  from  fracture  toughness  values  mea¬ 
sured  using  standardized  procedures. 

Both  Sumpter  [6]  and  Kirk  and  Dodds  [7]  achieved  good  agreement  between  the  initiation  frac¬ 
ture  toughness  of  single  edge  notched  bend,  SE(B),  specimens  and  structures  containing  part- 
through  semi-elliptical  surface  cracks  by  matching  thickness  and  crack  depth  between  specimen  and 
structure.  These  results  demonstrate  that  toughness  values  determined  from  shallow  cracked  SE(B) 
specimens  are  appropriate  for  assessing  the  fracture  integrity  of  structures.  However,  testing  stan¬ 
dards  have  not  yet  incorporated  formulas  permitting  evaluation  of  7  and  CTOD  for  shallow  cracks 
from  experimental  measurements  (i.e.  load,  crack  mouth  opening  displacement  (CMOD),  and  load 
line  displacement  (LLD)).  This  investigation  develops7  and  CTOD  estimation  procedures  applicable 
for  both  shallow  and  deep  crack  fracture  toughness  testing  for  materials  with  a  wide  range  of  strain 
hardening  characteristics. 


2.  APPROACH 

IWo  dimensional,  plane-strain  finite -element  analyses  of  SE(B)  specimens  are  performed  for 
crack  depths  from  0.05  to  0.70  a/H'with  Ramberg- Osgood  strain  hardening  coefficients  (n)  between 
4  and  50.  Tible  1  summarizes  the  conditions  considered.  The  analyses  provide  load,  CMOD,  and  LLD 
records  to  permit  evaluation  of  coefficients  relating/  and  CTOD  to  measurable  quantities.  The  range 
of  parameters  considered  in  these  analyses  allows  evaluation  of  the  dependence  of  these  coefficients 
on  alW  and  n.  The  estimation  formulas  divide  /  and  CTOD  into  small  scale  yielding  (SSY)  and  large 
scale  yielding  (LSY)  components.  In  each  formula,  the  SSY  component  is  defined  by  the  linear  elastic 
stress  intensity  factor,  JQ.  The  formulas  differ  only  in  the  L5y  component.  Procedures  to  estimate  the 
LSY  component  include: 

1 .  from  plastic  work  (area  under  the  load  vs.  LLDpi  curve,  or  Ap|  |  yj}) 

2.  CTOD/^  as  a  fraction  of  CMODpi  using  a  rotation  factor 

3.  CTOD/^.  from  plastic  work  (area  under  the  load  vs.  LLDp/  curve,  or  Ap|  |  yj,) 

4.  Jijy  and  CTODby  from  plastic  work  (area  under  the  load  vs.  CMODp/  curve,  or 
\iIcmod) 
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5.  CTOD/ijy  as  a  fraction  of  CMOU^/  without  the  notion  of  a  rotation  factor 
Existing  standards  employ  the  first  two  techniques  [1  -3];  the  remainder  are  new  proposals. 


Ihble  1:  SE(B)  specimens  modelled. 
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3.  J  AND  CTOD  ESTIMATION  PROCEDURES 

Current  Standards 

Existing  test  standards  for  J  and  CTOD  [1-3]  employ  the  following  estimation  formulas; 

K2(1-v2)  tip, 

J. — g —  +  _API 


pIlLLD 


K2(l  ~  v2)  r  )bCMOD., 


rp,b  +  a 


(3.1.1) 

(3.1.2) 


where 

K 

linear  elastic  stress  intensity  factor 

V 

Poisson’s  ratio 

"Ipi 

plastic  eta  factor 

B 

specimen  thickness 

b 

remaining  ligament,  W- a 

AplIlXD 

area  under  the  load  vs.  LLD^/  curve 

m 

constraint  factor 

°flow 

flow  stress,  average  of  yield  and  ultimate^ 

’■pi 

plastic  rotation  factor 

CMODp, 

plastic  component  of  CMOD 

Values  of  iip,,  m,  and  rp,  are  well  established  for  perfectly  plastic  materials  based  on  closed  form  solu¬ 
tions.  For  deeply  cracked  specimens  {a/W  ^  0.5),  current  test  standards  use  rip,  ■  2,  m  *■  2,  and 
rp,.)B0.44  Sumpter  [8]  and  Wu.,  et  al.  [9]  have  proposed  the  following  relations  to  account  for  crack 
depth  less  than  0.5  a/W: 

1.  ASTM  El  290  uid  BS  5762  both  uie  yield  itreu  in  the  CTOD  eitimation  equation.  In  thii  inveitiguion,  flow  itreii  is 
used  instead. 


2 


V-0.32+12^-49.5(^)' 

+  99.8(^)  for  a/W  <  0.282 

(3.1.3) 

Dpi  -  2.0 

for  a/W  0.282 

rp,  -  0.S  +  0.42^  -  4(^)“ 

for  a/W  <  0.172 

(3.1.4) 

tpi  0.463  —  0.04^ 

for  a/W  2  0.172 

Sumpter  derived  the  T)p|  equation  from  limit  analyses  of  the  SE(B),  while  Wu,  Cotterell,  and  Mai  used 
a  slip  line  field  analysis  to  determine  the  variation  of  rp|  with  alW.  Material  strain  hardening  alters  the 
deformation  characteristics  of  the  specimen,  thereby  altering  tip,,  m,  and  rp| .  Existing  procedures  ne¬ 
glect  any  influence  of  strain  hardening. 


3.2  New  Propoaale 

The  estimation  formulas  presented  in  Section  3.1  have  received  the  greatest  attention  as  the  coeffi¬ 
cients  relating 7  and  CTOD  to  experimental  measurements  are  amenable  to  closed  form  solution,  at 
least  in  the  non-hardening  limit.  For  hardening  materials,  closed  form  solution  is  not  possible,  there¬ 
fore  either  experimental  techniques  [10]  or  finite-element  analyses  [11]  are  used  to  provide  data 
from  which  t)p|,  m,  and  rp|  are  calculated.  Quantities  other  than  CMOE^/  and  Ap,  |  uj)  measured  dur¬ 
ing  a  test  can  also  be  related  to  7  or  CTOD,  if  the  proper  proportionaliQr  coefficient  is  known.  The 
following  are  some  alternatives: 


1. 


2. 


Estimate  CTODiQ,  from  plastic  work  (Apilnn)! 


_  K2(l  -  v2)  ,  tlc-L  . 


This  formula  is  analogous  to  eqn.  3.1.1  for  7  testing 


(3.2.1) 


Use  plastic  work  defmed  by  the  area  under  the  load  vs.  CMODn/  curve  (A,,.  I  rMnn)  to 
estimate  either 7i^  or  CTODfay:  ^  v-piicMOD/ 


K»(l-v2) 

E  ^  “Bb  Sil  CMOD 


(3.2.2) 


3. 


1 1 CMOD 


(.3.2.3) 


This  technique  eliminates  the  need  for  LLD  measurement,  which  simplifies  7  testing. 
Express  CTODi)^  as  a  fraction  of  CMOD^/: 

Eqn.  3.2.4  and  3.1.2  are  functionally  the  same,  thus  Tj^  and  rpj  are  related: 


Sorem  [11]  found  tpi  to  be  extremely  sensitive  to  the  CTOD -CMOD  relationship  for 
shallow  cracks.  This  estimation  procedure  was  proposed  to  circumvent  this  sensitivity. 
The  validity  of  this  approach  is  based  on  the  observed,  nearly  linear  dependence  of 
CTODbv  CMOD;,/  in  flnite -element  solutions. 
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In  this  investigation,  finite-clement  analyses  provide  data  from  which  Tip,,m,  rp,,  i1j_o  ’Ic-o 

and  Tig  are  calculated. 

4.  FINITE-ELEMENT  MODELLING 
TWo- dimensional,  plane  strain  finite-element  analyses  of  SE(B)  specimens  are  performed  using 
conventional  small  strain  theory.  The  analyses  are  conducted  using  the  POLO  -  FINITE  analysis  soft¬ 
ware  [12]  on  an  engineering  workstation. 

Uniaxial  stress  strain  behavior  is  described  using  the  Ramberg- Osgood  model 

where  Oo  is  the  reference  stress  (0.2%  offset  yield  stress  when  a  ■  1),  e©  *  Oo/E  is  the  reference 
strain,  o  ■  1,  and  n  is  the  strain  hardening  coefficient.  Strain  hardening  coefficients  of  4, 5, 10,  and 
50  model  materials  ranging  from  highly  strain  hardening  to  nearly  elastic  -  perfectly  plastic.  Figure 
1  illustrates  these  stress  -  strain  curves. 

J2  deformation  plasticity  theory  (nonlinear  elasticity)  describes  the  multi-axial  material  model. 
Tbtal  strains  and  stresses  are  related  by 


where  Sy  is  the  stress  deviator,  o*  is  the  Mises  equivalent  tensile  stress,  Oj^.  is  the  trace  of  the  stress 
tensor,  and  fly  is  the  Kronecker  delta. 


Figure  1:  Ramberg -Osgood  stress  strain  curves  used  in  the  finite -element  analysis. 


4 


Figure  2:  FlnitC’-element  model  of  the  a/(F>0.25  SE(B)  specimen. 

Finite- element  models  are  constructed  for  s/H' ratios  of  0.05,  0.15, 0.25, 0.50,  and  0.70.  The 
SE(B)  specimens  have  standard  proportions;  the  unsupported  span  is  four  times  the  specimen  width. 
Symmetry  of  both  geometry  and  loading  permit  use  of  a  half- symmetric  model.  Each  model  contains 
approximately  400  elements  and  1300  nodes;  the  ajW  ■  0.25  model  is  shown  in  Figure  2,  Eight- 
noded,  plane -strain  isoparametric  elements  are  used  throughout.  Reduced  (2  x  2)  Gaussian  integra¬ 
tion  is  used  to  eliminate  locking  of  the  elements  under  incompressible  plastic  deformation.  The  same 
half- circular  core  of  elements  surrounds  the  crack  tip  in  all  models.  This  core  consists  of  eight,  equal¬ 
ly  sized  wedges  (22.5  *  each)  of  elements  in  the  6  direction.  Each  wedge  contains  30  quadrilateral  ele¬ 
ments;  the  radial  dimension  decreases  geometrically  with  decreasing  element  distance  to  the  crack 
tip.  The  eight  crack-tip  elements  are  collapsed  into  wedges  with  the  initially  coincident  nodes  left 
unconstrained  to  permit  development  of  crack-tip  blunting  deformations.  The  side  nodes  of  these 
elements  are  retained  at  the  mid -point  position.  This  modelling  produces  a  Hr  strain  singularity  ap¬ 
propriate  in  the  limit  of  perfect  plasticity.  Crack-tip  element  sizes  range  from  0.2%  to  0.02%  of  the 
crack  length  depending  on  the  alW  modelled. 

Load  is  uniformly  distributed  over  two  small  elements  and  applied  at  the  center  of  the  compres¬ 
sion  face  of  the  specimen  to  eliminate  the  local  singularity  effects  caused  by  a  concentrated  nodal  load. 
Load  is  increased  in  30  to  50  variably  sized  steps  until  the  CTOD  reaches  5%  of  the  crack  length.  Strict 
convergence  criteria  at  each  step  insure  convergence  of  calculated  stresses  and  strains  to  the  third  sig¬ 
nificant  figure.  IV'o  to  three  full  Newton  iterations  at  each  load  step  are  required  to  satisfy  this  criteria. 
As  deformation  plasticity  is  strain  path  independent,  converged  solutions  are  load  step  size  invariant. 

The /-integral  is  computed  at  each  load  step  using  a  domain  integral  method  [13,14].  /  values 
calculated  over  domains  adjacent  to  and  remote  from  the  crack  tip  are  within  0.003%  of  each  other, 
as  expected  for  deformation  plasticity.  CTOD  is  computed  from  the  blunted  shape  of  the  crack  flanks 
using  the  ±  45*  intercept  procedure.  LLD  is  taken  as  the  relative  displacement  in  the  loading  direc¬ 
tion  of  a  node  on  the  symmetry  plane  located  approximately  0.46  ahead  of  the  crack  tip  and  of  a  node 
located  above  the  support.  This  procedure  eliminates  the  effect  of  spuriously  high  displacements  in 
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'nible  2:  Calculation  of  coefficients  in  J  and 
CrOD  estimation  formulas. 


Eqn. 

Coefficient 

3.1.1 

^Pi 

3.2.2 

ilj-c 

3.1.23.2.3 

3.2.13.2.4 

m 

3.2.1 

’IC-L 

3.2.3 

TlC-C 

3.1.2 

>-pl 

3.2.4 

% 

1 

•  j  (i  is  the  slope  of  this  line,  fp,  ■  l 

t:6-CTOD  ^ 


the  vicinity  of  both  the  load  and  support  points.  The  t),  m,  and  tpi  coefficients  are  determined  from 
these  results  by  calculating  the  slope  of  the  quantities  indicated  in  Ihble  2  at  each  load  step.  Slope  cal* 
culation  is  initiated  with  data  from  the  final  three  load  steps.  Data  from  earlier  load  steps  are  included 
in  this  calculation  until  the  linear  conelation  coefficient  (r)  falls  below  0.999.  This  procedure  elimi¬ 
nates  data  from  the  first  few  load  steps,  which  arc  predominantly  clastic,  and  therefore  not  expected 
to  provide  reliable  relationships  between  plastic  quantities. 


6.  RESULTS  AND  DISCUSSION 

The  variation  of  the  tj,  m,  and  rp|  coefficients  with  alW  and  n  determined  from  the  finite -element 
results  is  summarized  in  Figures  3— 4,  and  in  the  Appendix.  Solutions  for  non -hardening  materials, 
where  avaUable,  are  indicated  on  the  figures.  Each  coefficient  shows  considerable  variation  with  crack 
depth.  The  variation  with  material  strain  hardening  is  also  a  common  feature  of  all  coefficients  except 
hj-o  '''hich  relates to  Ap|  |  cmqd-  t1j-c  essentially  independent  of  «  for  a/W'a:  0.15.  The  remain¬ 
der  of  this  section  examines  the  differences  between  perfectly  plastic  and  finite -element  solutions, 
and  the  errors  associated  with  each  estimation  procedure.  Finally,  recommendations  of  /  and  CTOD 
estimation  formulas  for  use  in  fracture  testing  of  SE(B)  specimens  are  made. 
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5.1  Perfectly  Plastic  and  Finite  Element  Proportionality  Coefficients 
The  variation  of  both  tpi  and  T)j.with  alWiox  a  low  strain  hardening  material  (Figure  3  e-f)  agrees 
well  with  the  slip  line  field  solution  of  Wu,  et  al.  [9]  above  n/lf'ssO.lS.  However,  at  smaller  alW  the 
elastically  dominated  response,  ignored  in  the  slip  line  field  solution,  causes  a  deviation  between  the 
slip  line  field  and  finite-element  ipj  and  values. 

The  variation  of  T)p|  with  determined  by  finite-element  analysis  has  a  different  functional  form 
than  determined  by  Sumpter  [8]  using  a  limit  load  solution  (Figure  3a) .  The  limit  load  derivation  em* 
ploys  the  following  approximation  for  plastic  work: 


UpL  ■  PuM  ■  LLDpL 


(5.1.1) 


where 


PuM 


S 


5  ■  1  -  0.33^  - 
S  ■■  unsupported  bend  span 

Thus,  the  accuracy  of  t)p|  values  determined  by  limit  analysis  depends  on  the  equivalence  of  plastic 
work  calculated  by  eqn.  5.1.1  and  the  actual  plastic  work  (area  under  a  load  vs.  LLD^,/  diagram)  for 
a  strain  hardening  material.  This  equivalence  is  not  achieved  even  for  the  low  strain  hardening  n  >50 
material,  as  illustrated  in  Figure  5. 


5.2  J  and  CTOD  Estimation  Errors 

Figure  6  illustrates  the  variation  of/  and  CTOD  with  LLD  and  CMOD  for  an  a/if'">0.1S,  SH(B) 
determined  by  finite -element  analysis.  This  dependence  of  fracture  parameters  on  measurable 
quantities  is  contrasted  with  that  predicted  by  the  /  and  CTOD  estimation  procedures  using  t)  and  m 
coefficients  calculated  from  the  finite-clement  results.  Work-based  /  and  CTOD  estimates  (eqns. 
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0.000  0.005  0.010  0.015  0.020  0.025 

Plastic  Load  Line  Displacement  [inches] 

Figure  5:  Compariion  of  limit  solution  and  finite-element  reiulti  for  a/tt'aO.lS,  n  *50. 

3.1.1, 3.2.1,  3.2.2,  and  3.2.3)  match  the  finite-element  results  much  more  closely  than  do  formulas 
that  calculate  CroPi,^  as  a  fraction  of  CMOD^  (eqns.  3.1.2  and  3.2.4).  Figure  7  shows  /  and  CTOD 
estimation  errors,  more  clearly  illustrating  the  differences  between  the  estimation  procedures,  lb 
evaluate  the  effects  of  both  a/W  Mnd  n  on  estimation  accuracy,  the  following  error  measure  is  defined: 

ERR  »  ^  |E,FP*»,  /  (5.2.1) 

i-1  I  i-1 

where 

FpMt  _  ppft 

E,  ■ - - 1 100  percent  enor  at  load  step  i 

1  FP™, 

AT  total  number  of  load  steps 

FF***]  estimated  J  or  CTOD  at  toad  step  i 

FP^*i  /  or  CTOD  at  load  step  i  from  finite  -  element  analysis 

For  ano/ir*0.15  ln^5  SE(B),  theERRvalue  for  CTOD  estimation  using  r^i,  eqn.  3.1.2,  is  21%.  Com¬ 
parison  of  this  value  with  the  data  in  Figure  7demonstrates  that  ERR  is  a  root  mean  square  error  mea¬ 
sure. 

The  variation  of  ERR  with  aJW  and  n  for  the  six  estimation  procedures  is  shown  in  Figure  8.  Er¬ 
rors  associated  with  work-based  J  and  CTOD  estimates  (work  calculated  from  CMOD)  are  below 
5%  for  all  alW  and  n.  If  work  is  instead  calculated  from  LLD,  J  and  CTOD  estimation  errors  are  also 
generally  below  5%,  with  the  exception  of  shallow  cracks  in  a  very  low  strain  hardening  material 
(a/iFM0.05,  n  ibSO).  However,  equations  that  express  CTODiQ,  as  a  fraction  of  CMOE^/  are  inaccurate 
for  all  a/1V(EKR>n%)  in  highly  strain  hardening  materials  (n  s  S).  As  the  maximum  estimation  er- 
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2.0 


J  [In  •  klps/in^] 


0.C16 


CTOD  [inches] 


H  0.072  h 


0.008  h 


0.004  k 


0.000 


CTOD;,y  as  a  Fraction  ol  CMODp/ 
Eqns.  3.1.2  and  3.2.4 


Rnha  Elamtnt 
Results 


0.000  0.008  0.016  0.024  0.032 

CMOD  [inches] 


Load  Line  Diepiacement  [inches] 

Figure  6:  Viriation  of  7  and  CTOD  with  LLD  and  CMOD  for  a/IV-O.lS,  .n»5  SE(B). 


Figure  7:  J  end  CTOD  estimation  errors  for  n/W'-O.lS,  n»5  SE(B). 

ror  can  exceed  SJFOf  by  up  to  a  factor  of  2  (Figure  7),  EJiM  values  above  17%  are  clearly  excessive. 
Accuracy  improves  {EJlk<12%)  for  materials  with  less  strain  hardening  (n  a:  10).  However,  these 
estimates  have  accuracy  comparable  to  work-based  CTOD  estimates  only  for  deep  cracks  in  essen- 
tially  non-hardening  materials.  Thus,  the  validity  of  assumptions  made  in  deriving  the  various  es- 


Figure  8:  Variation  of  J  and  CTOD  estinnation  errora  with  a/W  and  n.  Symbols  represent  the  same 
conditions  in  each  figure,  (a)  eqn.  3.1.1,  (b)  eqn.  3.2.2,  (c)  eqn.  3.2.1,  (d)  eqn.  3.2.3,  (e)  eqn. 
3.1.2,  (f)  eqn.  3.2,4. 
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timation  procedures  directly  affects  their  accuracy,  J  and  CTOD  estimation  from  plastic  work  is 
achieved  by  partitioning  total  work  into  SSY  and  LSY  components,  Additive  separation  is  exact  be¬ 
cause,  for  a  linear  elastic  body,  K*(l  -  v)/Eistheclasticstrain  energy.  Conversely,  the  linear  relation 
between  CTODi^  and  CMOD^/  assumed  in  eqns.  3.1.2  and  3.2.4  cannot  exist  (exactly)  for  any  body 
with  an  elastic  component  that  varies  with  load  (i.e.  for  any  amount  of  strain  hardening).  Strain  hard¬ 
ening  strongly  influences  the  linearity  of  the  CTODi^  ->■  CMOE^/  relationship,  as  illustrated  in  Figure 
9.  Thus,  eqns.  3.1.2  and  3.2.4  work  best  for  minimally  strain  harde.ning  materials. 


CMODp/  [inches] 

Figure  9:  Effect  of  strain  hardening  on  the  linearity  of  the  CTOD/»  -  CMODw 
relation  for  a/WmO,SO, 


6.3  Recommended  J  and  CTOD  Estimation  Procedures 
6,3.1  Requirements  for  Accurate  Estimation 

The  formulas  used  to  evaluate  fracture  parameters  from  experimental  data  should  not  introduce  sub¬ 
stantial  errors  into  the  J  and  CTOD  estimates.  This  need  for  accuracy  favors  estimating  4^  and 
CTOD^  from  plastic  work.  Even  though  estimation  of  the  LSY  component  from  plastic  work  requires 
numerical  integ.  ation  of  experimental  data,  this  seems  warranted  to  reduce  errors  by  up  to  five~fold 
(compare  Figure  8d  to  Figure  8f).  In  addition  to  using  inherently  accurate  formulas,  selecting  t],  m, 
and  Tpi  coefficients  corresponding  to  a  specific  atW  and  material  should  not  be  a  potential  error 
source.  In  view  of  the  ambiguity  attendant  to  fitting  experimental  stress -strain  data  with  a  power  law 
curve,  insensitivity  of  t^,  m,  and  rp|  to  material  strain  hardening  would  be  extremely  advantageous. 
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5.3,2  J  Estimation 

The  only  procedure  that  meets  both  of  the  aforementioned  requirements  is  J  estimation  from  plastic 
work  based  on  CMOD.  By  fitting  the  data  in  Figure  3b,  the  variation  of  with  a/W  is  expressed 
as  follows: 

t),.c  ■  3.785  -  3.101^  +  2.018(^)*  for  all  n,  0.05  s  ^  s  0.70  (5.3.2.1) 

Figure  10  shows  this  fit  together  with  the  tij.^data.  The  use  of  values  from  eqn.  5.3.2.1  produces 

estimation  enors  of  at  most  9%,  and  generally  much  less,  as  illustrated  in  Figure  11.  In  situations 
where  fracture  toughness  in  terms  of  a  critical  J  value  is  desired,  estimation  using  eqns.  3.2.2  and 
5.3.2.1  is  clearly  superior  to  estimating /from  plastic  work  based  on  LLD,  where  qp|  depends  on  mate¬ 
rial  strain  hardening  coefficient.  Further,  estimating  /  from  CMOD  rather  than  LLD  eliminates  the 
need  to  measure  LLD,  which  simplifies  the  test  procedure. 

Despite  the  clear  advantages  of  estimating  /  from  plastic  work  based  on  CMOD,  estimation 
based  on  LLD  may  be  necessary  for  very  shallow  cracks  due  to  experimental  complexities  associated 
with  clip  gage  attachment  [15].  If/ estimation  using  LLD  is  unavoidable,  T)p|  can  be  indexed  less  ambig¬ 
uously  to  the  ratio  of  the  ultimate  strength  to  the  yield  strength  than  to  the  strain  hardening  coefficient. 
The  ultimate  tensile  strength  for  a  Rumberg-Osgood  material  is  obtained  by  solving  for  the  tensile 
instability  point,  converting  true  stress  to  engineering  stress,  and  taking  the  ratio  of  this  value  with 
0.2%  offset  yield  stress.  This  calculation  gives: 

The  variation  of  1/n  with  R  calculated  from  eqn.  5.3.2.2  is  shown  in  Figure  12.  This  figure,  along  with 
the  information  in  Tkble  Al,  is  used  to  determine  the  appropriate  T)p|  value  for  the  experimental  condi¬ 
tions  of  interest  based  on  data  from  a  simple  tensile  test. 


a/W 

'  Figure  10;  Comparison  of  eqn.  5.3.2.1  to  finite-element  data. 
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Figuri  11:  EiTor  aiiociated  with  uiing  t))  _  valuet  from  eqn.  5.3.2.I. 


Figure  12:  Relationihip  between  strain  hardening  coefficient  (n)  and  ultimate  to  yield 
ratio  (R)  for  a  Kamberg-Osgood  material. 


5.3.3  CTOD  Estimation 

As  noted  previously,  CTOD  estimation  from  plastic  work  is  considerably  more  accurate  than  CTOD 
estimation  directly  from  CMOE^/.  Use  of  eqn.  3.2.1  or  3.2.3  is  therefore  preferred  to  eqn.  3.1.2  or 
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3.2.4.  However,  the  ti,  m,  and  fpi  coefficients  in  all  of  these  equations  depend  strongly  on  n.  The  strain 
hardening  coefficient  is  estimated  from  R  as  described  in  section  5.3.1.  Appropriate  m  and  r|c_L  or 
T^c.cvalues  for  the  experimental  conditions  of  interest  are  then  determined  from  Ihbles  A3,  A4,  and 
A5,  respectively. 


6.  SUMMARY  AND  CONCLUSIONS 

Results  from  two-dimensional,  plane  strain  flnite-element  analyses  are  used  to  develop  /  and 
CTOD  estimation  strategies  appropriate  for  application  in  both  shallow  and  deep  crack  SE(B)  speci* 
mens.  Crack  depth  to  specimen  width  (a/JV)  ratios  between  0.05  and  0.70  are  modelled  using  Ram* 
berg-Osgood  strain  hardening  exponents  (n)  between  4  and  50.  The  estimation  formulas  divide  J  and 
CTOD  into  small  scale  yielding  (SSY)  and  large  scale  yielding  (LSY)  components.  For  each  case,  the 
SSY  component  it  determined  by  the  linear  elastic  stress  Intensity  factor,  A/.  The  formulas  differ  in 
evaluation  of  the  LSF  component.  The  techniques  considered  Include:  estimating  /  or  CTOD  from 
plastic  work  based  on  load  line  displacement  (Ap|  |  from  plastic  work  based  on  crack  mouth  open* 
ing  displacement  ( Ap)  |  cmod)*  plastic  component  of  crack  mouth  opening  displacement 

(CMOC^/).  Api  I  cmod  provides  the  most  accurate /  estimation  possible.  The  flnite-element  results 
for  all  conditions  investigated  fall  within  9%  of  the  following  formula; 

J  .  +  Sg^ApilcMOD  1J-c  - 3'/85’  - 

The  insensitivity  of  y)  j  _  c  to  strain  hardening  permits/  estimation  for  any  material  with  equal  accuracy. 
Further,  estimating/ from  CMOD  rather  than  LLD  eliminates  the  need  to  measure  LLD,  thus  simpli* 
tying  the  test  procedure.  Alternate,  work  based  estimates  for  /  and  CTOD  have  equivalent  accuracy 
to  this  formula;  however  the  t)  coefficients  in  these  equations  depend  on  the  strain  hardening  coeffl 
cient.  CTOD  estimates  based  on  scalar  proportionality  of  CTODSb^  and  CMOD^;  are  highly  inaccu* 
rate,  especially  for  materials  with  considerable  strain  hardening,  where  errors  up  to  38%  occur. 
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APPENDIX 

SUMMARY  OF  COEFFICIENTS  FOR  J  AND  CTOD  ESTIMATION 


Ikble  Al:  Variation  of  ijpi  with  a/W and  n  for  /  estimation  by  eqn.  3.1.1. 

•/W 

Ramberg~Osgood  Strain  Hardening  Coefficient  (n) 

4 

5 

10 

50 

0.05 

0.670 

0.746 

0.901 

1.192 

0.15 

1.295 

1.393 

1.542 

1.687 

0.25 

1.639 

1.686 

1.763 

1.753 

0.50 

1.924 

1.930 

1.924 

1.927 

0.70 

2.109 

2.130 

2.086 

2.052 

Tkble  A2:  Variation  of  iji/.c  with  a/W and  n  for  J  estimation  by  eqn.  3.2.2. 

a/W 

Ramberg-Oagood  Strain  Hardening  Coefficient  (n) 

4 

5 

10 

50 

0.05 

3.848 

3.793 

3.482 

3.420 

0.15 

3.359 

3J85 

3.322 

3.376 

0.25 

3.152 

3.138 

3.130 

3.1.17 

0.50 

2.748 

2.749 

2.728 

2,723 

0.70 

2.613 

2.641 

2.595 

2.562 

Uible  A3:  Variation  of  m  with  a/W  and  n  for  CTOD  estimation  by  eqni.  3.1.2, 3.2.1, 3.2.3, 
and  3.2.4. 

a/W 

Ramberg-Oagood  Strain  Hardening  Coefficient  (n) 

4 

5 

10 

50 

1.908 

1.786 

1.496 

1.291 

0.15 

1.963 

1.863 

1.573 

1.423 

0.25 

2.036 

1.938 

1.648 

1.501 

0.50 

2.177 

2.047 

1.788 

1.687 

0.70 

2.200 

2.093 

1.932 

1.810 
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'ftblc  A4:  Variation  of  i^c-l  with  atW  and  n  for  CTOD  estimation  by  eqn.  3.2. 1 . 


Ramberg-Osgood  Strain  Hardening  Coefficient  (n) 


0.335 


0.640 


0.795 


0.885 


0.959 


5 

10 

50 

0.402 

0.611 

0.800 

0.743 

0.982 

1.245 

0.872 


0.944 


1.018 


11ibleA5:  > 

/ariation  of  tjc-c  with  alW and  n  for  CTOD  estimation  by  eqn.  3.2.3. 

Ramberg-Osgood  Strain  Hardening  Coefficient  (n) 

a/W 

4 

5 

10 

50 

0.05 

1.929 

2.043 

2.310 

2.701 

0.15 

1.659 

1.806 

2.115 

2.493 

0,25 

1.530 

1.624 

1.904 

2.112 

0.50 

1.263 

1.344 

1.525 

1.605 

0.70 

1.187 

1.262 

1.341 

1.412 

TbbleA6:  > 

/ariation  of  rpi  with  a/W  and  n  for  CTOD  estimation  by  eqn.  3.1.2. 

a/W 

Ramberg-Osgood  Strain  Hardening  Coefficient  (n) 

4 

5 

10 

50 

0.05 

0.045 

0.053 

0.089 

0.142 

0.15 

0.132 

0.171 

0.261 

0.404 

0.25 

0.207 

0.240 

0.352 

0.431 

0.50 

0.292 

0.343 

0.380 

0.426 

0.70 

0.333 

0.341 

0.395 

0.398 

Tkble  A7:  Variation  of  rjt  with  alW  and  n  for  CTOD  estimation  by  eqn.  3.2.4. 


Ramberg-Osgood  Strain  Hardening  Coefficient  (n) 


0.459 


0.427 


0.382 


0.226 


0.125 


5 

10 

50 

0.499 

0.627 

0.729 

0.492 

0.595 

0.695 

0.418 

0.512 

0.563 

0.255 

0.274 

0.299 

0.127 

0.145 

0.146 
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